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ABSTRACT 

We review recent numerical studies and the phenomenology of spatially synchro- 
nized collective states in many-body dynamical systems. These states exhibit thermo- 
dynamic noise superimposed on the collective, quasiperiodic order parameter evolution 
with typically one basic irrational frequency. We concentrate on the description of the 
global temporal properties in terms of second-order difference equations. 



PACS: 05.45. 4-b, 47.25. Mr, 87.10.+e 



- 1 - 



Complex dynamical systems exhibit various types of temporal and spatial evolution 
ranging from uniform to chaotic. The mechanisms leading to self-organization mani- 
fested in constant (fixed-point), limit cycle, persistent-structure, quasiperiodic and other 
"collective" temporal behaviors have been a subject of extensive investigations. Un- 
fortunately, only for the simplest families of one-dimensional cellular automata, signif- 
icant progress in the understanding of individual rules has been achieved. A more 
phenomenological approach is based on classifying possible collective behaviors based 
on the experience with few-parameter systems such as iterative maps in one or more 
variables. The usual mean- field type idea is to lump the evolution of the complex system 
in one or more "collective" order parameters, for instance the total magnetization of 
an Ising-spin automaton, or generally large-wavelength Fourier components of the spin 
configuration. This approach, as well as the phenomenological characterization of the 
spatio-temporal patterns observed, are at the heart of most classifications of cellular 
automata.^~^° 

Recently, arguments have been advanced^ suggesting that spatially uniform col- 
lective evolutions in synchronous-dynamics complex systems are generically limited to 
"trivial" temporally constant or period- two cyclic states (at large times). The idea has 
been that for higher-period temporal evolution, with integer periods 3,4,... in units 
of the discrete updating time steps of the automaton, the domain walls between re- 
gions "out of phase" will have generically nonzero velocity. Thus, randomly generated 
small droplet regions "out of phase" with the surrounding configuration, are expected to 
spread irreversibly, destroying the spatial coherence. Similar arguments were advanced 
for quasiperiodic evolution, with domain walls replaced by regions of phase gradients. 

Surprisingly, numerical evidence has been reported indicating the existence of var- 
ious "nontrivial" self-organized states in higher-dimensional discrete-spin cellular au- 
tomata. Originally, these collective states were observed in certain d = 4, 5 totalistic- 
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rule automata. ° This discovery was followed by extensive numerical studies ' which 
both enriched the family of automata models showing these new collective states, and 
also classified various phenomenological properties to be listed shortly. Moreover, simi- 
lar collective states were found recently ■^^"■'^^ in continuous- variable coupled- map lattices 
in d = 2 through 6. Collective states were also observed in Id coupled- map lattices; 
however, their spatial synchronization pattern and other properties differ from collective 
states in d > 1. 

The models studied in Refs. 10-19 showed a great variety of self-organized collec- 
tive states (see especially Ref. 17). The most intriguing type is characterized by the 
formation of a stable loop-like structure in the return map, i.e., in the plot of the order 
parameter m{t + 1) vs. m{t). Figure 1 illustrates such a plot for the d = 5 totalistic 
automaton rule-^°'-^^ denoted R^g, in which a site is updated to the value crj(t -|- 1) = 1 
if the neighborhood sum over the von Neumann zone (the site itself and its 2d nearest 
neighbors) is between 5 and 9 at time t, and updated to ai{t + 1) = otherwise. The 
order parameter is defined as 

m{t) = J2<^jit)/j2^, (1) 

j j 

where the sum is over all the lattice points j, at which the spins crj{t) = 1 or are 
located. In Figure 1, about 1000 data points are shown after a transient of 5000 time 
steps has been discarded, for a system of 12^ sites. The initial condition was a randomly 
half- filled (with values 1) lattice; and boundary conditions were periodic. The large-time 
average value, M, was estimated and subtracted out. This quantity can be formally 
defined as 

M= lim 

T— >00 
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Bifurcation diagrams were studied as well, for coupled-map lattices for which there 
is a control parameter. ■'^^"■'^^ Restricting our discussion to the d > 1 models, we list sev- 
eral phenomenological observations based largely on numerical studies. Firstly, although 
noisy integral-period states (for example, period three) were found, noisy quasiperiodic 
behavior attracted more interest. Here we focus on such quasiperiodic states, typically 
with quasiperiods close to integral values 3,4,.. ., measured, for example, by the winding 
numbers^^"-^-*^'-^^ about a point near the origin in plots like Figure 1. 

The noise exemplified by the plot in Figure 1, was found to decrease "thermody- 
namically" with increasing number of spins, according to 



noise amplitude ~ 1 




This important observation indicates that the spatial synchronization mechanism is 
statistical, short-range in nature, rather than global, i.e., that the synchronization is 
achieved as a local effect, via correlations of adjacent system parts. Unfortunately, the 
nature of the spatial self-organization has not been understood theoretically thus far. 

Regarding the temporal evolution, numerical evidence based on maximum-entropy 
analyses, as well as some phenomenological considerations, -"^^ suggest that the appar- 
ently quasiperiodic behavior observed in the new collective states both in d — 1 and 
d > 1, has in fact periodic time dependence with a single basic frequency. The quasiperi- 
odicity is thus solely due to the incommensurability of the associated period with the 
clock period At = 1 of the automaton time steps. This observation can be quantified 
by the relation 



m{t » 1) = M -h ttn sin ( + 0n j + noisy contributions . (4) 
n=l ^ ^ 



The noisy contributions have no characteristic frequencies associated with them, and 
they are presumably important only for d > 1 models, as mentioned earlier. 

All the above properties were obtained for fully deterministic evolution. However, 
stability of the collective states was checked^^'^^ by adding a small stochastic element 
in the evolution rules and for a range of initial conditions. Empirically, these states 
disappear when the stochastic element increases, or when the initial conditions signif- 
icantly depart from the random uniform configuration. In the latter case the system 
evolves into a trivial, typically dilute, disordered-type state. Configurations with im- 
posed domain walls, and the way the system "heals" such perturbations, were studied 
as well.^^ 

The theoretical situation remains largely unsatisfactory. Indeed, with the accu- 
mulated numerical information, it would be of interest to understand both the spatial 
and temporal self-organization mechanisms of these new collective states. One step in 
this direction, focusing on phenomenological conclusions one can draw from the numer- 
ical observation (4) , was the proposed second-order in time nature of the discrete-time 
dynamics of the collective behavior at hand.^^ In the rest of this Short Review we 
summarize the approach of Ref. 16. 

Empirically, one observes that the leading oscillatory term in (4) is usually the 
dominant contribution: 

|ai| > \an>i\ . (5) 

Such time evolution for t » 1 can be described by the linear second-order dynamical 
rule 



H{t) = 2 cos{c)n{t - 1) - n{t - 2) , 
- 5 - 



(6) 



where 



fM{t) = m{t) -M . (7) 

Indeed, the solution of (6) is 

//(t)=aisinf— + , with ^= — • (8) 

The linear evolution rule (6) only contains one parameter, c, determining the basic 
frequency. Linear fits of the plots of the discrete second-order derivative, //(t + 1) — 
2ji{t) + iJi{t — 1), vs. for several different models, -"^^ yielded period values consistent 
with both the maximum-entropy and winding-number estimates. Figure 2 illustrates 
such a plot for the data used in Figure 1. 

The noisy curve in Figure 2 is not precisely a straight-line, nor even single- valued. 
However, it is well known that such "Lissajous-like" distortions can be incorporated by 
allowing for nonlinear terms in the dynamical rule, which we now rewrite as follows. 



oo 

H{t+l)-2n{t) + n{t-l) = -4sin2(c/2)//(t) + ^6„//"(t) + noisy contributions. (9) 

n=2 

The form of the contributions leading to noise in finite-size samples is not clear at 
present. 

In some cases, the single-harmonic "linear" second-order difference equation (6) 
works quite well. As an illustration, we present in Figure 3 a plot analogous to Figure 2, 
but for a d = 3 automaton described in Ref. 13: the spin is updated to 1 if the sum 
over the 3d von Neumann zone is or 5, and set to otherwise. However, examples 
are known of large nonlinearity^^ as well as of "Lissajous-shaped" features obscured by 
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the noise, especially in high dimensions where the system sizes accessible to numerical 
simulations are small. 

A direct maximum-entropy fit of the parameters of the first six terms in (4), i.e., 
including the terms up to gq but excluding the noisy contribution, yields the n{t) se- 
quence which generates the continuous curve in Figure 4. If we only keep the leading 
harmonic, (8), then the elliptical-shaped curve is obtained; see Figure 4. The general 
equation of such an elliptical-shaped single-harmonic return map is 



+ 



1 2 



2 cos(c) 



= sin^(c) 



(10) 



While these observations shed little light on the nature of the spatial self-organization, 
they indicate that when the noise can be neglected, i.e., presumably in the infinite- 
system limit, the collective states essentially follow a "mechanistic" time evolution which 
in terms of the average order parameter m(t) is fully deterministic and reversible . The 
latter property follows from the second-order (in time) nature of the dynamics. As a 
result, the droplet- type arguments which are largely irreversible in nature, should be in- 
correct although the mechanism of their breakdown can be fully illuminated only after 
the spatial self-organization is described theoretically. Loosely speaking, reversibility 
implies that for any locally disordering mechanism there should exist a reversed locally 
ordering mechanism. It seems likely that neither of these processes involves droplets. 

In summary, we surveyed a new emerging field of study of collective, self-organized 
behavior in complex systems. The phenomenological information collected by extensive 
numerical simulations suggests that these systems offer an interesting challenge of theo- 
retical description of spatially self-organized and synchronized evolution, for which the 
conventional equilibrium "droplet" stability concepts do not apply. While the spatial 
mechanisms remain to be explained, the time-dependence at least of the global order 
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parameter, could be quantified phenomenologically by the second-order discrete-time 
dynamics outfined in this Short Review. 

The authors gratefully acknowledge discussions and collaboration with B. Buck and 
V.A. Macaulay. This research was supported in part by the Science and Engineering 
Research Council (UK). 
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FIGURE CAPTIONS 

Fig. 1: Return map, m{t + 1) — M vs. m{t) — M, for the cellular automaton rule R^g, 
described in the text: about 1000 time steps are shown. A transient of the first 
5000 time steps was discarded. System size was 12^. 

Fig. 2: Second-order discrete time derivative, the left-hand side of (9), plotted vs. //(t), 
for the data shown in Figure 1. 

Fig. 3: Second-order time derivative plot similar to Figure 2, for the 3-dimensional 
Hemmingsson model, described in the text. 

Fig. 4: The data of Figure 1 are shown together with their "noiseless skeleton" recon- 
struction obtained by keeping the first six oscillatory terms in (4), i.e., up to gq. 
The ellipse corresponds to keeping only the leading oscillatory term, corresponding 
to relations (6), (8), (10). 
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